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ABSTRACT 

Diffraction  radiation  is  that  electromagnetic  energy 
which  is  caused  by  a  relativistic  charged  particle  passing 
through  an  aperture  in  an  opaque  material.   Ter-Mikaelian 
solved  for  the  diffraction  radiation  from  a  point  charge. 
This  paper  discusses  the  phenomena  resulting  from  finite, 
relativistic  charge  bunches. 

Using  the  Huygens-Fresnel  principle,  diffraction 
patterns  from  spherical  and  cylindrical  charge  distributions 
are  found  and  plotted.   For  charge  bunch  sizes  less  than  the 
radiation  wavelength,  the  results  are  almost  identical  to 
those  for  point  charges. 

The  radiation  pattern  is  composed  of  two  regions.   The 
"transition  region"  is  characterized  by  a  strong  peak  at 
0  =  y""'"'  "^he  Lorentz  factor.   The  "diffraction  region" 
consists  of  a  series  of  peaks  and  nulls  in  field  strength 
typical  of  the  standard  plane  wave  diffraction  pattern. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION  7 

A.  BACKGROUND 7 

B.  THE  VIRTUAL  PHOTON 8 

C.  DIFFRACTION  THEORY  BASED  ON 

HUYGENS'  PRINCIPLE  10 

D.  FRESNEL  AND  FRAUNHOFER  DIFFRACTION  12 

II.  COMPUTATION  OF  DIFFRACTION  RADIATION  16 

A.  SOLUTION  FOR  A  POINT  CHARGE 16 

B.  SOLUTION  FOR  A  SPHERICAL  CHARGE  BUNCH  18 

C.  DIFFRACTION  RADIATION  FROM  A  LARGE 

SPHERE  OF  CHARGE 25 

D.  DIFFRACTION  RADIATION  FROM  A  CYLINDRICAL 

CHARGE  BUNCH  31 

E.  APPROXIMATIONS  AND  LIMITATIONS  38 

III.  PARAMETERS  AFFECTING  DIFFRACTION  RADIATION  39 

A.  DEPENDENCE  OF  DIFFRACTION  RADIATION  ON  Y  4  0 

B.  DEPENDENCE  OF  DIFFRACTION  RADIATION 

ON  BUNCH  RADIUS 43 

C.  DEPENDENCE  OF  DIFFRACTION  RADIATION  ON 
APERTURE  RADIUS  45 

IV.  DIFFRACTION  RADIATION  POWER  49 

V.  CONCLUSIONS 51 

APPENDIX  A:   EVALUATION  OF  LINE  CHARGE  APPROXIMATION 

FOR  CYLINDRICAL  BUNCH  53 

APPENDIX  B:   DIFFERENCES  WITH  THE  THEORY  OF 

TER-MIKAELIAN  58 

LIST  OF  REFERENCES 61 


BIBLIOGRAPHY  

62 

INITIAL  DISTRIBUTION  LIST  

53 


I.   INTRODUCTION 

A.   BACKGROUND 

Diffraction  radiation  occurs  when  a  relativistic  charged 
particle  passes  through  a  hole  in  a  conducive  or  dielectric 
material.   This  phenomenon  is  very  well  understood  for  single, 
infinitesimal  charges.   Ter-Mikaelian  computed  the  diffraction 
radiation  for  a  single,  relativistic  electron.   [Ref.  1]   Less 
understood  are  the  effects  of  charge  bunching.   Presumably, 
bunch  sizes  on  the  order  of  the  radiation  wavelength  will 
significantly  affect  the  diffraction  radiation  intensity  by 
destructive  interference.   This  paper  discusses  diffraction 
radiation  from  charge  bunches  with  radii  greater  than  zero  but 
less  than  the  radiation  wavelength.   Diffraction  radiation 
will  be  computed  for  spherical  and  cylindrical  bunch  shapes. 
In  addition,  the  effects  of  bunch  size,  beam  energy,  and 
aperture  radius  will  be  explored. 

When  electrons  are  travelling  at  relativistic  velocities, 
they  can  be  treated  as  "pseudo-photons."   A  method  very  simi- 
lar to  that  used  in  physical  optics  based  on  the  Huygens 
principle  is  then  used  to  find  the  diffraction  radiation  pat- 
tern.  [Ref.  l:p.  376]   Panofsky  and  Phillips  describe  the 
"virtual  photon  concept"  in  Reference  2. 


B.   THE  VIRTUAL  PHOTON 

The  electric  field  from  a  "fast  electron"  approaches 
that  of  a  plane  wave  as  the  electron's  velocity  approaches 
c,  the  speed  of  light  in  a  vacuum.   Representations  of  the 
E-fields  of  electrons  at  various  velocities  are  shown  below 


V=r  0 


V    <  c       v^c       v-C 


For  a  uniformly  moving  electron,  one  can  solve  the 
Lienard-Wiechert  potentials  to  find  the  field  as  viewed  by 
a  stationary  observer.   One  finds  the  component  of  the  elec- 
tric field  perpendicular  to  the  electron's  motion  is  very  much 
greater  than  the  parallel  component  and  is  given  by 


Ej_{t) 


y 


eb 


4tt£   '  ,k2  ^  /  ^  ^2.  3/2 
o   (b  +  (vty)  ) 


where : 


e 
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=   the  electron  charge; 

=   the  perpendicular  distance  from  the 
electron's  path; 


the  Lorentz  factor  =  ( 
v/c . 


^/^T^^)-^ 
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One  can  view  this  field  as  a  sharp  pulse  of  electromagnetic 
energy.   The  magnitude  of  this  pulse  is  highly  dependent  on 
Yf  which  defines  the  relativistic  energy  of  the  electron. 
In  addition,  the  duration  of  the  pulse  decreases  as  the 
electron  energy  rises.   The  electric  field  of  an  untra- 
relativistic  electron  will  appear  very  nearly  as  a  delta 
function  in  the  time  domain.   A  representation  of  this 
electromagnetic  pulse  is  shown  below: 


The  Fourier  transform  of  E, (t)  gives  the  frequency  components 
of  this  pulse. 


^Goi 
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where 


Integrating,  Panofsky  and  Phillips  found 


E  .    =    —^ — n^)K,  A] 

o 

where  K,  is  the  modified  Bessel  function  of  the  second  kind. 
K, (x)  approaches  infinity  as  1/x  near  x  =  0  and  goes  to  zero 

_  y 

as  e    for  arguments  greater  than  one.   Panofsky  and  Phillips, 
as  well  as  Ter-Mikaelian  in  his  development  of  diffraction 
radiation,  approximated  E  .  by  equating 


for      (jo    < 


_  ;    417    ebv 

OJl 

£  .  yv 

for      CO    <    ~-    . 
b 


In  his  treatment  of  diffraction  radiation  from  a  single 
particle,  Ter-Mikaelian  assumes  that  the  Fourier  components 

of  E  .  are  spatially  limited  to  a  circle  of  radius  b  z    Ay. 

col       '^       -^  ' 

He  states  that  the  particle  will  radiate  at  a  frequency  co 

only  when  Ay  >  r  ,  where  r   is  the  characteristic  dimension 
^  '     o         o 

of  the  aperture.   [Ref.  l:p.  377] 

C.   DIFFRACTION  THEORY  BASED  ON  HUYGENS '  PRINCIPLE 

Ter-Mikaelian  used  Huygens '  principle  to  compute  the 
diffraction  radiation  pattern  from  a  charged  particle  passing 
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through  an  opening.   He  thought  of  each  point  in  the  aperture 

as  a  source  of  secondary  waves  and  then  integrated  over  the 

aperture  to  find  the  total  field  at  some  distant  observation 

point.   The  field  at  each  point  within  the  aperture  is 

assumed  to  be  unaffected  by  the  obstruction.   In  practice, 

edge  effects  will  be  apparent  for  small  apertures. 

Huygens  Principle  states  that: 

every  point  on  a  primary  wavefront  serves  as  the  source 
of  spherical  secondary  wavelets  such  that  the  primary 
wavefront  at  some  later  time  is  in  the  envelope  of  these 
wavelets.   Moreover,  the  wavelets  advance  with  a  speed 
and  frequency  equal  to  that  of  the  primary  wave  at  each 
point  in  space.   [Ref.  3:pp.  60-61] 

Using  this  theory,  one  can  describe  the  propagation  of 

a  wave  past  an  obstruction  or  aperture  by  mathematically 

summing  the  contributions  of  each  "Huygens  wavelet."   Fresnel 

expanded  upon  Huygens'  principle  by  stating: 

every  unobstructed  point  of  a  wavefront,  at  a  given  instant 
in  time,  serves  as  a  source  of  spherical  secondary  wavelets 
(of  the  same  frequency  as  the  primary  wave) .   The  ampli- 
tude of  the  optical  field  at  any  point  beyond  is  the 
superposition  of  all  of  these  wavelets  (considering  their 
amplitudes  and  relative  phases).   [Ref.  3:p.  330] 

There  are  some  problems  in  using  the  Huygens-Fresnel 

principle,  most  notably  the  effects  of  electron  oscillators 

at  the  edge  of  the  aperture.   Remember  that  Ter-Mi]<:aelian 

assumed  that  the  field  within  the  aperture  was  unaffected  by 

the  obstruction.   Oscillating  electron  "clouds"  at  the  edge 

of  the  aperture  will  produce  their  own  contributions  to  the 

total  electric  field  within  the  aperture.   The  field  is  then 

made  up  of  a  part  from  the  incident  electric  wave  and  one  from 
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the  induced  oscillations  at  the  edge  of  the  obstruction. 
However,  for  apertures  larger  than  the  wavelength  of  the 
incident  radiation,  Huygens-Fresnel  diffraction  theory  is 
expected  to  provide  quite  accurate  results. 

D.   FRESNEL  AND  FRAUNHOFER  DIFFRACTION 

If  the  observation  plane  is  located  relatively  far  from 
the  aperture,  one  can  ignore  the  vector  aspects  of  wave 
propagation  and  use  the  scalar  Huygens-Fresnel  principles. 
For  each  small  source  at  the  primary  wavefront,  the  contribu- 
tion to  the  field  at  observation  point  P(X,Y)  is  given  by 


dE  {X,Y)   =   E  (x,y)  exp  i(wt-kr)  ds/r 
D  a 


The  exponential  term  provides  the  phase  information  and 
the  r  term  in  the  denominator  takes  into  account  spherical 
spreading.   A  diagram  of  the  problem  is  shown  below: 


AferturC 
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In  the  typical  diffraction  problem  encountered  in  the 
study  of  physical  optics,  a  plane  wave  is  assumed  at  the 
aperture,  so  the  electric  field  across  the  aperture  has  no 
spatial  variation  in  amplitude,  polarization  or  phase.   When 
dealing  with  the  radiation  from  charged  particles,  however, 
one  must  consider  that  the  field  is  radial  with  amplitude 
changing  with  distance  from  the  charge.   In  our  problem,  we 
locate  the  observation  point  on  the  X-axis  and  deal  only 
with  the  x-component  of  the  field  at  the  aperture.   The  field 
is  symmetric  about  the  Z-axis. 

The  total  field  at  P(X,0)  is  given  by  the  superposition 
of  all  the  Huygens  wavelets  from  all  the  sources  at  P(x,y) 
of  size  ds  across  the  aperture. 


r.   /w  n\         r     1- /    N  sxo  i(ojt-kr)  -, 
E   (X,0)   =      j      E(x,y)  — ^ ds 

^  t^^  Aoerture 


This  is  a  very  difficult  integral  to  evaluate  in  that  the 
exponential  term  has  a  modulus  of  27t  ,  making  it  difficult  to 
make  an  asymptotic  approximation. 

The  phase  factor  can  be  approximated  by  taking  the  binomial 
expansion  of  r  and  approximating  it  to  the  first  order. 


99  9  999  22 

r    =   Z   +  (X-x)   +  (0-y)    =   Z^  +  X""  -  2Xx  +  x^  +  y 


9     2  2     2 

Z^  +  X  )  -  2Xx  +  x^  +  y 


R^(l  -  2Xx/R^  +  (x^  +  y^)/R^) 
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where 


2       2     2 
R    =   X   +  Z 


so 


2      2     2    2  1/2 
R(l  -  2XX/R   +  (x   +  y^)/R'^)^/^ 


The  binomial  expansion  is  given  by 


/I  ,   \  n      T  ,     ,  n  (n-1)   2  , 
(1  +  x)    =   1  +  nx  +  — ^^-j — -   X   + 


To  first  order. 


2     2 

r   =   R  -  Xx/R  +  (x   +  y  )/2R 


The  total  field  at  P(X,0)  becomes 


E    (X,0)   =   e^(a3t-KR)     J     ^a^x^^-y^ 
^  Aperture      r 


^  ^iK((x^+y^)/2R  -Xx/R)  ^g 


For  large  R,  one  can  replace  the  r  in  the  denominator  with 
R  and  the  integral  becomes 

l(wt-KR)  -ik(— -^-— ) 

Ep^^(X,0)   =   ^ ^ /E^^(x,y)e  ^ds. 
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If  R  is  much  greater  than  the  aperture  size  and  the 

phase  difference  from  all  of  the  incremental  sources,  ds , 

2   2 
across  the  aperture  is  negligible,  then  the  k(x  +y  )/2R 

term  in  the  exponential  is  very  small  and  may  be  neglected 

In  general,  if  this  term  is  less  than  one  radian,  then  the 

integral  can  be  approximated  by 


i  (ojt-KR) 

E^  ^(X,0)   = /e   (x,v)exp(ikXx/R)ds 

p ,  X  K      ^  ap{ 


In  the  Fraunhofer  approximation  the  electric  field  at  a 
distant  point  is  just  the  two-dimensional  Fourier  transform 
of  the  electric  field  in  the  aperture. 
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II .   COMPUTATION  OF  DIFFRACTION  RADIATION 

A.   SOLUTION  FOR  A  POINT  CHARGE 

A  single  relativistic  electron  has  a  Fourier  spectrum, 
as  derived  earlier,  given  by 

e    r  /ojb^  ^.  ,a)b 


E,,   =   2^ [(^)K,  (^)J 


"wl      .    2      ,  vv   1  VY 

4tt  r  bv    '       ' 


o 


Using  the  Fraunhofer  method,  the  field  at  a  distant  observa- 
tion point  is  given  by 


^  Aperture 


where  the  x-component  of  the  field  at  the  aperture  is  given 
by 


E^  X   =  -2^ [(^)K,  (^)j  cos  e  . 


Changing  to  polar  coordinates  gives  us 


i(ajt-KR)      O     27T  ,        ,  -iv^     n/T. 

e  f  f  e        r/wbsT^    /wbM  ^      ikXbcose/P.  ,     ,„  -,, 

z /        /      TT— -tW^^lW^^"^^®  ^  ^^^"^ 

0         0         4tt  e  bv       '  ' 

o 

i  (cot-KR)  O       ,  ,  27T 


/      Ak,  (^)      /       cosee^^^^/^)^°^^   d6db 

-*  VY       1    VY       J 


4/e  V  0  VY     1  VY     Q 

o 


16 


Integrating  over  dQ  yields 


i  (cot-KR)  o 


e  e       _     _a3_       j       bK,  (^^  ^    '^^ 

^     J  1  vy 


47T£    V  "^        n  1    VY        1      R 


i  (ojt-KR)        .  /I  (jo/1 

^  le  ,     CO   V     r  O  JV      ^T        /I     .       ^\Tr        /  0\ 

Z 2?F^   V  ^TV— TIITT^  R>^2(^or)^i(V 


yv     1  ^      o  R     2     yv 


which  contains  the  final  form  of  the  diffraction  radiation 
pattern.   As  Ter-Mikaelian  did  in  Reference  1,  the  radiation 
field  is  found  by  subtracting  the  particle  field  from  the 
results  above.   The  particle  field  is  found  by  letting  t 
go  to  infinity  in  the  integral  over  the  aperture.   It  is 
interesting  to  note  that  as  i   goes  to  infinity,  the  diffrac- 
tion radiation  field  strength  goes  to  zero.   Physically,  at 
this  point  the  electron  no  longer  "feels"  the  presence  of  the 
obstruction. 

Ter-Mikaelian,  in  Reference  1,  computes  the  diffraction 
radiation  from  a  charged  particle  passing  through  a  circular 
aperture.   Using  the  method  outlined  in  the  previous  section, 

he  determines  (for  ai   <<  1) 

o 


E  ^(x,o)  =  ^ o-^-T  J  (q-^^) 

P,X  27T£C   2     2   o  ^  o 
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where 


0) 

a   =   —  ; 
vy 

q   =   k  sin  0  =   the  projection  of  the  wave  vector 

k  on  the  plane  z  =  0 ; 

n,        =      the  aperture  radius. 
o  ^ 


This  result  appears  to  be  significantly  different  from  the 
expression  derived  earlier.   In  reality,  however,  the  two 
plots  are  virtually  identical.   (See  Appendix  B.) 

B.   SOLUTION  FOR  A  SPHERICAL  CHARGE  BUNCH 

The  diffraction  problem  has  been  solved  for  single 
charged  particle  radiation. 

The  process  for  a  finite  charge  bunch  can  be  solved,  as 
a  first  approximation,  by  assuming  a  spherical  charge  dis- 
tribution of  radius  a   with  constant  charge  density,  p. 

In  the  frame  of  the  moving  charge,  the  electric  field 
is  radial  with  a  magnitude  given  by 


1 

r  p   for   r  <  a 


E   = 


3£_  ^  ^   ^^^   ^  -  "o 


^o  P    , 

for   r  >  a 


2j- 1^  J.    J.  "       a. 
—  o 
J£  r 
o 


This  result  is  easily  obtained  using  Gauss'  law.   A  repre- 
sentation of  E(r)  is  shown  below. 


To  find  the  field  from  the  moving  charge  distribution  in 
the  frame  of  a  stationary  observer,  one  must  find  the  Lorenz 
transformation  of  the  radial  field.   With  a  spherical  charge 
distribution,  an  observer  views  a  field  emanating  from  the 
center  of  the  charge  bunch;  the  field  from  the  uniform 
spherical  distribution  appears  to  be  a  point  charge  equal  to 
the  total  charge  of  the  sphere.   If  the  observer  is  located 
within  the  sphere,  he  still  "sees"  a  point  source,  but  of  a 
reduced  strength  (that  portion  of  the  sphere's  charge  enclosed 
by  a  Gaussian  surface).   Consequently,  the  Lorenz  transfor- 
mation of  the  moving  spherical  charge  bunch  must  have  a  form 
similar  to  that  of  the  point  charge.   In  the  lab  frame,  one 
finds  a  perpendicular  electric  field  of 


Ej_  (t) 


4/3  TTa  p 
o 

47Te_ 


4/3nb  p 
4tt£ 


"^b 

T — ^ — T — T^      tor   b  >  a„ 

,  ,  ^  ,   ,  -2,  2 ,  3/2  -   o 

I  [vt.)      +  Y   b  ]  ^ 


-^ 7^ — =i — :^-7-^      for   b  <  a 

2  .  -2,  2 ,  3/2  -  o 


[(vt)"  +Y  "b"] 
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This  form  of  E  (t)  is  only  approximate  for  b  -^^  a  .   For 
small  bunches  (a   <<  1  ) ,  the  error  introduced  by  this 
approximation  is  negligible. 

Since  E  (t)  of  the  sphere  has  the  same  time  dependence 
as  that  of  the  point  charge,  it  must  have  Fourier  components 
of  the  same  form. 


Uw 


f  4/3TTa    p         ,  , 

.    2      ^    ^Vy^^IVy^J 
4tt    £    bv         '  ' 

o 


4/37Tb    p  r  ,ojb,  ,,     |^^^^ 


k   47T     £     bV 

o 


for      b    >    a 
—      o 


for      b    <    a 
—      o 


Again,  E  .  is  negligible  for  A  <  b/Y ,  we  make  the  additional 
requirement  that  the  wavelength  be  large  compared  to  the 
bunch  size  in  order  that  the  bunch  radiate  coherently. 

Again,  one  proceeds  using  the  Huygens-Fresnel  principle 
to  find  the  radiation  pattern  from  the  spherical  bunch  passing 
through  a  circular  aperture.   With  the  observation  point 
located  on  the  X-axis,  one  expects  to  see  only  an  X-component 
to  the  perpendicular  field.   Due  to  symmetry,  only  the  x- 
component  of  the  field  across  the  aperture  contributes  to 
the  total  field  at  P(X,0).   The  transverse  electric  field 
at  P(X,0)  is  given  by 


i  (cot-KR) 


E    (X,0,Z) 
P  /  A 


/ 

Aperture 


E^^^(x,y)e^^^/^  dS 
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where    the   x-component   of    the    field   at    the    aperture    is    given 
by 


3 

4/3TTa    p         ,  , 

— T [( — )  K-,  ( — )]cos    9      for      b    <    a 

4^2,    ^^      vy       1    vy 

E  ^  ° 

a,X 


4/3TTb    P^(a3b)K,  (^)]coi 

^      2         ,  VY         1     VY 


.  -  ,  .       .  . .       for      b    <    a 

.      -         ,  VY         i     VY  —        o 

4tt    £    bv  ' 


Changing  to  polar  coordinates  yields 


i  (ojt-KR)  .  ,,  ..,      „  ,^> 

E_  ,,   =   ^ . /E    rb,e)e^(^^^^°^  Q/^^  d: 


p,X  Z      ^     a,X 


where 


ds   =   b  dedb 


Finally,  one  obtains 


^Ko^t-KR)  (b,e)ei^^^^°^Q/^^bdedb 

Z  -*    a,X 


i  (cot-KR)  o       2Tr    .  /-  ,3 


^____     r  r       4/3Trb  p     ojb         a3b     ^^^  g^KkXbcos  e/R)^  dedb 

2  0^        0^        4Abv     ^     1^ 

o 


0         4^2^  bv    ^     1  ^ 
o  o 
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We  shall  first  solve  for  the  integral  from  b  =  0  to  a  ; 

o 


we   have 


o       2tt    .  ,_  ,3 


f         J       VStt^       a)b         o^b     ^^^     ^i  (kXbcos  0/R)   ^  ^ 
0  0         4Abv         ^     1^ 


a  9 

O  TT 

•^  b    ( — )Kt  ( — )  cos  ee  ^       dSdb   . 

^  VY       1    V7       J 


/I     2  -^  VY       1    VY        „ 

4Tr  £  V  0  0 

o 


The  integral  over  d0  is  equal  to  (2TTi)J-,  (kXb/R)  .   The  integral 
becomes 


a 
o 


^^.    4/3Trp  .  ,  3,a)b,^,     /Ojb^^     ,kXb.      ,, 

2ui         -^  /  b     (--)K,  (— -)  J.  (-5-)     db    . 

.2  ^J                     Vy        1     VY        1        R 

4-rT    £„v  0                       '               ' 


'O 


As  an  example,  one  can  expect  to  see  parameters  in  an 
experiment  on  diffraction  radiation  in  the  micro-wave  region 
as  follows: 

f      =      8    GHz 

a        =      0 . 5    cm 
o 

Y     •         =20 
'mm 


v      ~      c=3xl0         cm/s 


k      =      27T/A      =      1.675    cm    ^ 


X/R   =   sin  9  where  9  is  the  angle  between  k  and  the 
Z-axis . 
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Using  these  values,  the  maximum  values  of  the  arguments  in 
K-|  (x)  and  J-,  (x)  in  the  integral  from  0  to  a   are 

(wb/yv)   =   8.33  x 10~^ 

(kXb/R)   =   8.38  X 10"^ 

With  these  values  for  arguments,  the  two  Bessel  functions  in 
the  integral  can  be  very  closely  approximated  by 


K^(x)   =   1/x 


and 


J-^(x)   =   x/2 


The  integral  becomes 


a  5 

T    .    4/3TTP      f       ,3  1  ,kXb.  -,,  „    .      4/37TP  .1  „.  ,  o. 

2^^  Wi^    I       ^    2  ^^^^     =     2TT1  -f-^(2  K)  (-^)    sin 
o      0  4tt    £  V 

o 


We  next  evaluate  the  integral  from  a   to  t  .   Since  the 

^  o      o 

integral  here  has  the  same  e-dependence  as  in  the  previous 
case,  the  integral  becomes 


4/37Ta    p  o        ,  ,  1  .., 

2TT1   —2 /  (y^)Ki(— )Ji(-^)    db    . 

47T    e    V      a  '  ' 

o  o 
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Note  that  in  this  case,  we  have  not  made  the  approximations 
for  the  two  Bessel  functions  since  b  will  in  practice  extend 
to  values  much  larger  than  one  wavelength.   This  integral 
can  be  easily  solved  and  yields 


o 


/      W^^lW^^l^^^^     ^      Skx,2      ,03,2^ 


r^    r  ,kX>-,    ,     o,  ^    ,       o>       w    ^    ,       o...    -     o. 
^    [^o^  (-R ^ ^1  ( V "^2  (-R-)   -^  '^l  (-^) ^2  ( V 


,  ,,         Loa  kXa  kXa  o^ 

r,KX>^^    ,     Ov  T-    -        o,  w    T    ,        o,,.,    .     o,-,. 


Combining  the  solutions  of  the  two  integrals  and  sub- 
tracting the  particle  field,  the  total  transverse  field  at 
P(X,0)  is  given  by 


3  2 

l(cot-KR)  4/3713    P  T     w       ^  1 

^p,x(^'0)      =     ^-^— (2.i^-^){{|(Jf)(-|)  .^(  ^/  ,) 

47T    £V  (^)       +  (— ) 


,^         oj/L  kXl  kXY  L0I. 


,  „         coa  kXa  kXa  oja 


k  sin  9  1 


(ksinP)^  +{—)^ 
VY 
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C.   DIFFRACTION  RADIATION  FROM  A  LARGE  SPHERE  OF  CHARGE 

For  diffraction  problems  where  the  radiation  contribution 

from  within  the  charge  bunch  is  significant  (i.e.,  a   ~  A  ), 

^  ^  o     o 

one  can  not  make  the  approximation  used  in  the  previous  section 

where  E,  had  the  same  time  dependence  for  b  <  a   as  it  did  for 
i  ""^  o 

b  >  a  .   In  this  section,  we  shall  gain  a  more  exact  under- 
o  ^ 

standing  of  the  problem  by  determining  the  diffraction  radia- 
tion from  a  large  sphere  of  charge. 
Using  Gauss's  law,  one  finds 


I 


p o  b 

4 11  e      2   r 


for   r  >  a 
—   o 


p   4/3TTr   b 
2    r 


4tt£ 


for   r  <  a 
—  o 


where  b/r  provides  the  perpendicular  component  of  E,  as  shown 
below: 


■zr^ 
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This  relation  reduces  to 


P       <  ' 


4/3TTa^b 

T —        for      r    >    a 

3  o 


4tt  £ 


4/3   Tib  for      r    <    a 


o 


A   representation    for    E      is    shown   below    for    several    values    of 
b: 


b^  a. 


/1^<>t*«/«a+»»'»  for    b^C\„ 


Figure    1 


This  function  for  b  <  a   can  be  integrated  numerically 
only.   However,  one  can  make  a  very  close  approximation  for 
E   with  b  <  a   by  writing 
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4/37Tb(1.5a  -Z) 


F 


4TTe 


1.5a 
:  o 

<     4/3^b 


n/272 
-  Va  -b 
o 


4/3TTb(1.5a  +Z) 


1.5a 


-  Va  -b 


a  -b 
o 


for 


for 


for 


v^a^-b^   <   Z   <   1.5a 
o  o 


-  V  a  -b     <   Z   <     Va"-b 
o  o 


-1.5a       <   Z 
o 


x/  2  ,2 

va  -b 

o 


o 


This  approximation  is  shown  on  Figure  1.   For  most  values 

3 
of  b,  this  linear  approximation  to  the  1/r   dependence  of 

E| (r)  outside  the  sphere  lies  very  close  to  the  actual  curve 

This  is  particularly  M:rue  as  b  approaches  a  ,  where  most  of 

the  diffraction  radiation  from  within  the  sphere  is  expected 

to  originate.   Replacing  z  with  -vt ,  E  (t)  is  equal  to 


(  4/3Tra  b 
'        o 


Ej_(t) 


(b2.(vt)2)3/2 


4/3^b(1.5a  +vt) 


4tt£ 


1.5a 
o 

4/37Tb 


y  2 , 2 

-  V  a  -b 
o 


4/3^b(1.5a  -vt) 


1.5a 


-  Va^-b^ 


a  -b 
o 


for     b  >  a 


r~2     2 
for     b  <  a   ,  va  -b     <  -vt  <   1.5a 
o         o  o 


A~2     2  Hi     2 

for     b<a,-va-b     <  -vt  <  v  a  -b 
o  o  o 


12     2 
for     b  <  a   ,   -1 . 5a     <  -vt  <  -  v  a  -b 
o  o  o 


For  a  relativistic  electron  bunch,  the  Lorenz  transforma- 
tion gives 


27 


f    4/37Ta  b 
o 


E    (t)      = 


4Tr£ 


Y    < 


,^2    ,    ^,2.3/2 
(b  +(YVt)    )    ' 


4/3Trb(1.5a  +Yvt) 


-  va  -b 


1.5a     -va-b 
o  o 


4/3'rTb 


4/3TTb(1.5a  -Yvt) 


I,    1 .  5a     -    va-b 
o  o 


-    va  -h 


for     b  >  a 


o 


y  2  .2 

va  -I 


a^-b"  1.5a 

for     b  <   a   , '-  -  s/ 1:  <  ° 


o'  Y 


Y 


7  2,2  J  2  /I 

for     b  <  a   ,  ^ <  -vt  < 

o'  Y  Y 


J  2  ,2 
-va-b 

for     b  <  a   ,  °  <  -vt  <  ° 


-1.5a 


o'        Y 


Y 


or 


E    (t)      = 


r  4/3^a  b 


(b^+(Yvt)^)^/^ 


4/3TTb(1.5a  +Yvt) 


v/  2  ,  2 
-     va-b 


4tt£ 


o 


1.5a 
o 

4/3ub 


4/3TTb(1.5a  -YVt) 


1.5a 


o 


■-  v/?:? 


o 


for     b  >  a 


o 


for     b  <  a   , 


-1.5a 


o'         vy 


for     b  <  a   ,     t     < 
o      '    ' 


y  2 , 2 

°<   t<-  ° 


VY 


va  -b 
o 

VY 


7  2,2 
va-b 

for  b  <  a   , <  t  < 

o  V  Y 


1.5a 


VY 


The  diffraction  radiation  found  by  integrating  over  the 
aperture  froir.  a   to  t   for  the  large  sphere  of  charge  will 
be  identical  to  that  found  in  the  small  sphere  approximation 
and  need  not  be  discussed  here. 
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This  section  will  discuss  only  that  radiation  originating 
from  the  aperture  for  0  <  b  <  a  .   For  this  range  of  b,  the 
Fourier  components  of  the  field  are  found  by  writing  E.  (t) 
as  a  difference  of  two  simple  functions  whose  Fourier  trans- 
forms are  tabulated  in  any  introductory  text  on  Fourier 
analysis.   E. (t)  is  written  as  the  difference  of  two  tri- 
angular pulses,  as  shown  below: 


It  is  easily  seen  that 

^     1.5a 
(     4/3TTb(1.5  a^)        sin    (ar.^— ^) 


I     1    c  ,/  2  ,  2  0) 

y     1 .  5a    -  V  a  -h 


o  o 


4/37Tb(1.5a   ) 
-  4( )    -  4/3TTb) 


sm 


2  My 


y  2 , 2 

va  -b 
o 


1.5a    -  ^a^-b^ 
o  o 


(""-2W)     \ 


0) 


r 
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A  diagram  of  a  typical  contribution  to  E    from  each  function 
is  shown  below: 


(jj 


<,lr 

H'->' 

2  V^ 

SM 

L 

One  can  easily  see  that  E^^  ^  is  cut  off  at  high  frequencies 

The  field  at  a  distant  observation  site  P(X,0,Z)  is 
given  in  the  usual  manner  by 


E  y(X,0,Z)   = 

p,  A 


i  (Wt-KR) 


27T 


/    /   E   cos  ee 


KkXbcos  e/R) 


0    0 


coi 


b  dbdG 


i(ojt-KR) 


2TTi   / 

0 


o 


\.'i'^'^^- 


This  integral  can  be  evaluated  numerically  to  find  the  dif- 
fraction radiation  pattern  from  a  large  sphere  of  charge. 
When  this  is  done,  the  results  are  not  accurate  enough  to 
make  detailed  predictions  regarding  the  diffraction  pattern, 
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They  do  show,  however,  that  the  diffraction  pattern  from  a 
large  sphere  of  charge  is  consistent  with  that  found  for 
small  bunches.   The  field  strengths  are  significantly  less, 
but  the  pattern  of  peaks  and  nulls  is  roughly  the  same. 

D.   DIFFRACTION  RADIATION  FROM  A  CYLINDRICAL  CHARGE  BUNCH 

Consider  first  a  stationary  charge  distribution  with 
axes  as  shown  below: 


>    ?(^''''^) 


We  wish  to  find  the  electric  field  at  an  observation  point 
P(X,Y).   By  cylindrical  symmetry,  the  electric  field  at  an 
observation  point  on  the  X-axis  will  have  only  i  and  k  com- 
ponents.  We  are  interested  primarily  in  the  i  component 
since,  for  the  relativistic  electron  bunch,  Ej_  >>  E..  . 

In  order  to  solve  for  Ei ,  one  must  first  find  the  electro- 
static solution  for  a  cylindrical  bunch.   One  can  make  an 
approximation  by  treating  the  cylindrical  bunch  as  a  line 
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charge  of  linear  charge  density  p  (ira  ).   For  b  <  a  ,  the 

o  o 

2 

linear  charge  density  is  given  by  p  (irb  )  .   A  diagram  of  the 

problem  is  shown  below: 


For   b    >    a 


o 


E(X,Z) 


pTia 


?/ 


4^e    '    (x2    4-(Z-Z')^^/^ 


dz' 


Let    u   =    Z-z;    du    =    dz . 


pTia 


o 


'^^    Z-L        (b24u2)3/2 


du     = 


u 


b      (b  +u  )  Z-L 


pira 


o  1 


Z-L 


""'^   "^    (b2+z2)V2        ,^2^,2-L)2)V2' 


For    a    slow  moving   bunch,    located    at    the    aperture    at    t    =    0, 
one   can    substitute     (-vt)     for    Z    and   write    the    time   varying 
field    at    the   aperture    as 


pTTa 


Ej_{t) 


o 


(-vt) 


vt+L 


4TT£b    ,,2,    .,2.1/2        ,,  2,  ,   .  ,_  ,2,1/2 
(b  +(vt)    )    ^  (b  +(vt+L)    )   ^ 
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For  a  relativistic  charge  bunch,  the  component  of  the 
electric  field  perpendicular  to  the  bunch  path  is  given  by 


E  (t)  = 


I 


4Treb  ^  \,  2_^  .  ^  .2.1/2  ^  ,,2^,^  ^T^2^1/2   ^^'^  ^   ^o 
(b  +(vtY)  )      (b  +(vtY+L)  )  ^ 


P  o  ^  ,     -vty      ,    (vty+L)       ^    , 

47TE:b    ,,  2  .  ,  ,2.1/2   ,,  2 ,  ,  .  ,^,2,1/2  o 

(b  +(vtY)  )  ^    (b  +  (vty+L)  )  ^ 


As  in  the  previous  section,  the  Fourier  components  of 
the  field  for  b  >  a   are  given  by 


oo 

— OO 


DTTs               *°                    icot                             °°           icot 
1    ^     o      r      r      (vty+L)e ,.  r     vtye -,.  . 


Evaluating  the  first  integral,  one  finds 


°°    ,   .     ,^.    iojt  °°    , vty+L.    loot 

\4)^+ (vty+L)  ^  ^        ^/^    ,  ,vty+L  " 

"■  b 


letting  ijj  =  (vty+L) /b;  t  =  (bi|;/vy )  -  (L/vy )  ;  and  di|j  =  (vy/b)dt, 
the  integral  becomes 
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°°  iaj( — ^b .coL       1 — ip 

b    r  ,     VY    VY  V.   1 —    *^  I   VY 


—00 


/ 2  "^      -oo    / 2 


di|; 


■j  Q 

Using  the  equation  e    =  cos  6  +  i  sin  9  and  remembering  the 
properties  of  even  and  odd  Fourier  integrals,  the  integral 
is  rewritten  as 


o,   1 —   "^  i[isin( — t^)        .2b   vy  ^ .  -oob. 


Using  the  same  techniques  to  evaluate  the  second  integral 
in  E   ,  one  finds 

CO  J. 


2  .ojL 

i  ^  ^  (2b)  J,    (a3b)  (^V  _i)        for     b  >  a 

2tt  4Tr£b  '    vy     1  vy  —    o 


E 

00 


"L 


„  .ooL 

7^—  7 — r-Y  ( — )K,  ( — )  (e         -1)        for     b  <  a 

27T  47T£b       vy     1  vy  —    o 


As  previously  explained,  the  Hankel  transform  of  the  field 
at  the  circular  aperture  is  evaluated  to  find  the  field  at 
observation  point  P(X,0). 
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i  (wt-KR) 
^  e '_  _i_ 

■'p,X  Z  2tt 


O  ZTT  2        ^,  ,  1 

f  r        pTTb  ^,  ,2b,  ^    ,wb,  ,     VY   -, ,         „ 

^i        J       4iIb^(^)K^(^)(^     ^-DcosOe 


ikXbcost 


'^^^bdbd( 


^  ^  2 

O        27T    pTTa 


•  coL 
1- 


+       /         /       __gY(^)K,(^)(e^-l)coseei^^^°^Q/^^bdbd^ 


a  0 

o 


47T£b     vy     1  vy 


Simplifying,  the  equation  becomes 


i(ajt-KR)   .        ^    i— 


(  f       ,3^    ,cob,    ikXbcos  0/R  ^.j, 

b  K,  ( — )e  ^      cos  ec3b 


de 


9 

+  a^   /  °  /  \k,  (!^)ei^^^^^^  e/I^  cos  edb  de 
o 


As  in  the  problem  for  a  spherical  charge  bunch,  integrating 
over  d9  gives 


^i(ojt-KR)       _   i — 


-^ 


f       ,3„  /tob,  ^  ,kXbvj,  ,2   r   UT^  /'^\ -r  /J'^^\  ji. 
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Evaluating  the  integral  from  0  to  a  ,  one  makes  the  usual 
approximation 


a 

o 


since  in  this  range,  K  (x)  1/x.   One  finds 


a  P     kXa 

Yi    /  °  b2j^(3s|b)^^    ^    ^te^  ^o  ^2(-R^: 

CO  ^^        1   R  o) 


Next,  evaluate  the  integral  from  a_  to  '^  . 

^  o  o 


A. 
2         f        ,^^    /Wb.  ^     ,kXb>  -,,  2  1 

^o         /        ^^1^^^^1^-T-^^^      =      ^ 


^    VY       1      R  o     ,kX>  2    .    ,   0)  .  2 

^O  ^-r)       ^^^ 


,  ^         kX^  (jj^  kX-l  03-^ 


1  ..        kXa  coa  kXa  coa 

-    ^o^-T^2(-R-^^1^^     -    ^  ^l^-R-)^2(^^^ 


Since  X/R  =  sin  9,  the  total  field  at  the  observation  plane 
is 
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■  /   j_      T^\  2      .coL 

1  (ojt-KR)    DTTa        1 — 

E    (X,0)      =     ^—- l-^(e^^-l) 

p  Z  27Te;v 


^         —  V  "^-    Q  Jo  (k  sin  e  a  ) 
I  I     GO   K  Sine     2  o 


( 


(  (k  sine)^  +{~) 


U  {k  sinej„(k  sine.'L  )K^  (— ^)    -  —  J^  (ksinB^  )K„(— ^)} 
o  2  olvy  vyl  o2vy 


oja  coa 

-a  {k  sinej^(k  sinGa  )K,  ( — -)    -  —  J,  (k  sinGa  )K„( — -)])] 
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The  effect  of  bunch  length  is  in  the  (e   ^-1)  term.   Expand- 
ing the  exponential  in  this  term  gives  an  imaginary  and 
real  part  of  the  solution. 


.a3L          i(^)/2     i(^)/2    -i(^)/2 
(e  ^^-1)   =   e   ^^    2i  ^ —2^ 


i(^)/2 
->  •     VY       .   ,  coL  V 
2ie         si^^I^) 
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At  first  glance,  this  function  seems  to  give  the  diffraction 
an  oscillating  dependence  on  radiation  frequency  for  all 
frequencies.   One  must  remember,  however,  that  the  K,  (coa  /vy ) 
term  in  the  solution  very  nearly  cuts  off  the  radiation  at 
to  >  vy/a  .   The  real  and  imaginary  terms  in  this  expression 
give  an  additional  phase  factor  to  the  overall  solution  for 
diffraction  radiation.   For  wavelengths  much  larger  than 
bunch  size,  the  real  term  is  approximately  zero  and  the 
dependence  on  bunch  length  is  given  by  i ( — ) . 

E.   APPROXIMATIONS  AND  LIMITATIONS 

The  reader  is  cautioned  to  remember  that  in  making  the 
Fraunhofer  approximation  to  the  diffraction  problem,  one 
assumes : 

1)  The  diffraction  angle  is  small  (<  1  rad) . 

2)  The  phase  difference  from  all  "sources,"  ds ,  across 
the  aperture  is  small  (k(x^  +y2)/2R)  <  1). 

3)  The  aperture  size  is  greater  than  the  wavelength  of 
the  incident  radiation.   (Neglect  edge  effects.) 

4)  Since  charge  was  assumed  to  be  uniformly  distributed 
throughout  the  bunch,  the  quantized  nature  of  charge 
was  ignored.   This  can  be  done  if  the  electron  den- 
sity is  sufficiently  high  (n"-'-/-^  <  X). 

5)  Charge  bunch  sizes  are  assumed  to  be  less  than  the 
radiation  wavelength  to  ensure  coherent  radiation  of 
the  bunch.   Larger  bunches  have  smaller  intensities 
due  to  destructive  interference. 
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III.   PARAMETERS  AFFECTING  DIFFRACTION  RADIATION 

In  the  previous  chapter,  diffraction  radiation  was 
computed  and,  for  a  given  frequency,  was  found  to  be  depen- 
dent upon  several  parameters,  namely  y ,    bunch  radius,  ao , 
and  aperture  radius,  ^  .   In  addition  to  these  effects,  the 
radiation  is  significantly  cut-off  at  high  frequencies. 

Remembering  that  both  the  small  sphere  and  cylindrical 

charge  bunches  had  Fourier  components  proportional  to 

(cob/vy)  KCcob/vy),  one  can  make  the  approximation  that  the 

fields  are  small  for  w  >  vy/b.   (When  the  argument  of  K-,  (x) 

is  1,  xK(x)  =  .6019.   This  gives  a  field  intensity  of  about 

one  third  the  maximum  value.)   If  the  maximum  value  of  b  is 

the  bunch  radius,  a  ,  then  the  maximum  value  of  oj  is  given 

o  ^ 

by  CO     =  vy/a  .   One  expects  radiation  at  2aj     to  be  10  dB 
-'   max     "  o  ^  max 

less  than  the  peak  values. 

Figure  2  can  be  used  as  a  rough  guide  to  the  frequency 
dependence  of  observed  radiation. 


Figure  2 
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For  the  large  sphere,  the  Fourier  components  of  the 
fields  were  found  to  drop  off  as 


sm  ( .  75 


VY 


2 

0) 


One  can  say  that  the  field  is  cut-off  at  co  =  vy/a  . 

It  is  easily  seen  that,  for  each  bunch  shape  discussed, 
the  diffraction  radiation  is  very  near  the  peak  values  at 
microwave  frequencies.   A  previous  experiment  at  the  NPS 
Linac  discovered  diffraction  radiation  in  the  microwave 
region  (f  =  8  GHz)  from  a  circular  aperture  [Ref.  4].   For 
f  =  8  GHz,  cob/vyK-i  (ajb/vy)  is  very  nearly  equal  to  unity. 

The  plots  of  diffraction  radiation  vs.  9,  the  off -axis 
angle,  consist  of  two  distinct  regions.   The  first  is  charac- 
terized by  a  strong  peak  at  9  =  y   •   Since  the  peak  at 
9  =  Y    is  typical  of  transition  radiation,  this  region  of 
the  curve  is  called  the  "transition  region." 

At  9  >>  Y   /  the  radiation  pattern  appears  more  as  the 
typical  diffraction  pattern  of  peaks  and  nulls.   This  region 
is  called  the  "diffraction  region." 

A.   DEPENDENCE  OF  DIFFRACTION  RADIATION  ON  y 

The  y-dependent  portion  of  the  function  describing  dif- 
fraction radiation  tends  to  be  very  complicated  and  hard  to 
understand  analytically.   The  strongest  dependence  is  in  the 
transition  region,  where  a  very  strong  peak  in  field  strength 
occurs  at  9  =  y   .   In  this  region,  the  field  strength  is 
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2  7 

very  nearly  given  by  (k  sinO)/((k  sinO)   +(w/vy)  )).   The 

strength  of  this  peak  is  directly  proportional  to  y-   In 
the  diffraction  region,  the  effect  is  less  pronounced. 
One  expects  to  see  more  diffraction  radiation  energy  from 
electrons  of  higher  energy.   Figure  4  shows  the  dependence 
of  diffraction  radiation  on  y.   The  curves  for  a  cylindrical 
bunch  will  be  identical  to  this  figure.   Increasing  y  tends 
to  sharpen  the  diffraction  pattern,  increasing  the  strength 
of  the  peaks  in  field  strength. 

For  pseudo-photons  of  a  given  frequency,  the  transverse 
fields  tend  to  extend  farther  out  for  increasing  y.   Alterna- 
tively, the  field  for  a  given  frequency  and  radial  distance, 
b,  tends  to  be  higher  for  an  increasing  y.   Figure  3  demon- 
strates this  point.   If  one  compares  the  same  Fourier  com- 
ponents of  two  relativistic  electrons  with  different  energies, 
the  more  energetic  electron's  fields  extend  much  further  out 
than  those  of  the  less  energetic  electron. 


^^^) 


Figure  3 
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Ter-Mikaelian    determined    that,    for    a   point    particle 
with    WL   /vY    <<    1/    the    ydependence   was    given   by 


k    sin 


(k  sin0)^   +  (— )^ 

VY 


The  Y-dependence  found  in  the  last  chapter  for  both  cylin- 
drical and  spherical  charge  bunches  is  consistent  with  that 
of  Ter-Mikaelian.   One  expects  the  fields  to  go  to  zero  for 
non-relativistic  particles. 

B.   DEPENDENCE  OF  DIFFRACTION  RADIATION  ON  BUNCH  RADIUS 

Throughout  this  paper,  we  have  assumed  that  the  bunch 
radius  is  smaller  than  the  wavelength  of  the  radiation.   When 
this  is  the  case,  one  may  assume  coherent  radiation  from 
the  bunch. 

By  examining  the  equations  for  diffraction  radiation  for 
all  three  bunch  types,  one  can  see  that  diffraction  radiation 
tends  to  decrease  for  bunches  with  larger  radii.   For  small 
bunch  radii,  a  general  decrease  in  the  field  strengths 
from  the  peak  levels  for  infinitesimal  bunch  sizes  is  observed, 
as  shown  in  Figure  5.   For  cylindrical  bunches  of  equal 
charge,  the  diffraction  radiation  field  strength  is  inversely 
proportional  to  L,  where  L  is  the  bunch  length.   This  is 
expected  since  the  radial  fields  at  the  aperture  are  propor- 
tional to  the  charge/length  within  the  bunch.   The  (e      -1) 
factor  also  specifies  a  length  at  which  coherent  radiation 
from  the  bunch  can  be  expected. 
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One  can  easily  understand  the  effect  of  bunch  size  by 
examining  the  static  radial  fields  from  a  sphere  of  charge. 
The  effect  from  a  cylindrical  charge  bunch  will  be  similar. 


Figure  6 

Figure  6  compares  the  fields  from  two  spheres  of  radii 

a   and  a'.   As  shown  earlier,  for  an  infinitesimal  sphere, 

the  fields  increase  assymptotically  to  infinity  as  r  goes 

to  zero.   For  a  sphere  of  radius  a  ,  the  field  increases 

^  o 

with  the  same  1/r  dependence  as  r  decreases  from  infinity 
to  ao  and  then  drops  off  linearly  to  0  at  r  =  0 .   One  can 
easily  see  from  the  figure  that  the  sphere  of  smaller  radius 
is  the  source  of  more  electrostatic  energy  than  the  larger 
sphere.   As  the  sphere  gets  smaller,  the  total  electrostatic 
energy  rapidly  approaches  that  of  a  point  charge. 

C.   DEPENDENCE  OF  DIFFRACTION  RADIATION  ON  APERTURE  PJ^DIUS 
A  plane  wave  incident  upon  a  circular  aperture  results 
in  the  familiar  Airy  diffraction  pattern.   The  large  central 
maximum  in  intensity  is  known  as  the  Airy  disk.   The 
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diffraction  radiation  pattern  from  a  relativistic  electron 
beam  is  very  similar,  except  that  the  central  maximum  is 
replaced  by  a  null.   This  central  null  results  from  the 
radial  nature  of  the  electron  fields.   In  the  Fraunhofer 
diffraction  method,  integrating  over  the  aperture  for  dG 
results  in  a  0th  order  Bessel  function  of  the  first  kind 
for  the  plane  wave  and  a  1st  order  Bessel  function  for  the 
electron  diffraction  radiation,  thus  explaining  the  princi- 
pal difference  in  the  two  radiation  patterns. 

Figure  7  shows  the  diffraction  radiation  pattern  for 
cylindrical  and  small  spherical  charge  bunches  passing 
through  circular  apertures  of  various  radii.   Changing  the 
aperture  size  affects  the  diffraction  radiation  pattern  in 
both  the  transition  and  diffraction  regions.   One  can  define 

a  dimensionless  parameter,  k-^  ,  where  k  =  — k  =  the  wave 

'^  o  c 

vector  and  n.      is  the  aperture  radius.   This  parameter  alone 

determines  the  location  of  the  peaks  in  the  diffraction  region 

When  k't   increases,  the  peaks  gain  in  strength  and  move  in 

closer  to  the  0=0  axis.   Decreasing  k^   spreads  the  dif- 

^    o   ^ 

fraction  pattern  out  over  wider  angles.   This  is  predicted 

since  the  diffraction  radiation  pattern  is  caused  principally 

by  the  J,  ^  (k-1  sin  9)  terms  in  the  solution.   Note  that 
^  1 , 2    o 

patterns  with  a  large  number  of  maxima  and  minima  are  found 
only  for  aperture  sizes  much  larger  than  the  wavelength  of 
the  observed  radiation.   This  is  to  be  expected  since,  for 
the  normal  diffraction  of  a  plane  wave,  a  small  aperture  acts 
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almost  as  an  isotropic  point  source  of  radiation.   In  the 

case  of  diffraction  radiation,  we  still  see  the  null  at 

0=0  due  to  the  radial  nature  of  the  fields.   Since  the 

dimensionless  parameter  is  given  by  ki  ,  one  can  change  the 

appearance  of  the  diffraction  pattern  by  changing  .%    ,    the 

physical  size  of  the  aperture,  of  k,  which  determines  the 

wavelength  of  the  radiation. 

In  the  transition  region,  changing  the  parameter,  k^    , 

affects  the  strength  of  the  peak  at  0  =  y   •   As  ki   goes 

to  zero  the  peak  at  y    increases,  asymptotically  approaching 

the  value  characteristic  of  transition  radiation.   As  k^i 

o 

increases,  the  strength  of  this  peak  decreases,  until,  as 
k^      goes  to  infinity,  the  peak  goes  to  zero.   At  this  point, 
the  charge  no  longer  senses  the  obstruction  and  no  energy  is 
radiated . 

One  should  also  notice  that  the  diffraction  radiation 
pattern  for  spherical  and  cylindrical  charge  bunches  are 
identical  for  bunch  lengths  less  than  the  radiation  wave- 
length.  For  very  large  charge  bunches,  interference  effects 
are  expected  to  reduce  the  overall  intensity  of  radiation/ 
(bunch  charge) ,  but  the  location  of  maxima  and  minima  will 
not  change.   The  locations  of  the  radiation  pattern's  maxima 
and  minima  are  principally  a  function  of  radiation  wave- 
length and  aperture  size,  as  it  is  in  the  normal  plane  wave 
diffraction  problem. 
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IV.   DIFFRACTION  RADIATION  POWER 

It  has  been  suggested  that  diffraction  radiation  may- 
be used  for  electron  beam  targetting  or  in  a  non-interceptive 
beam  diagnostic  device.   [Ref.  5] 

Present  diagnostic  systems  are  based  on  the  beam  passing 
through  a  thin  foil  and  using  transition  or  Cerenkov  radia- 
tion to  analyze  beam  characteristics.   The  problem  with  this 
scheme  is  that  the  beam  is  disturbed  by  the  foil.   If  the 
beam  passes  through  an  aperture,  one  can,  in  theory,  obtain 
the  same  diagnostic  information  without  disturbing  the  beam. 
If  diffraction  radiation  is  to  be  used  as  a  diagnostic  or 
targetting  means,  one  must  gain  an  appreciation  for  the  power 
levels  expected  from  a  typical  device. 

For  an  electromagnetic  wave  given  by 


E   =   g^e^^^t-^^)  , 


the  power  per  unit  area,  or  irradiance,  is  given  by 


CE    o 

<S>   E   ^  E^ 
2    o 


One  can  compute  the  diffraction  radiation  irradiance  by 
using  this  expression  together  with  the  equations  derived 
in  Chapter  II. 
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As  an  example,  the  diffraction  radiation  power 
delivered  by  a  1000  amp,  50  MeV  electron  beam  passing 
through  a  20  cm  diameter  circular  aperture  will  be  maximized 
at  6  =  10\   If  one  assumes  10  cm  long  cylindrical  bunches 
of  radius  .1  cm  and  evaluates  the  expression  for  diffraction 
radiation  from  a  cylindrical  charge  bunch,  the  irradiance 
at  Z  =  1  m  from  the  aperture  is  equal  to 


T   =  -^   E^   =   5  00  KW/m^  . 
2    o  ^ 


The  power  per  unit  solid  angle  is  500  KW/SR. 
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V.   CONCLUSIONS 

In  this  paper,  the  equations  regarding  diffra-:tion 
radiation  from  a  relativistic  charge  distribution  passing 
through  an  aperture  have  been  derived  from  the  Huygens- 
Fraunhofer  treatment  of  the  fields  across  the  aperture.   The 
calculations  are  very  similar  to  the  familiar  problem  of 
plane  wave  diffraction  solved  in  every  physical  optics 
class,  with  the  complications  that  the  fields  are  radial 
with  varying  field  strength  in  space. 

While  a  plane  wave  incident  on  a  circular  aperture  causes 
a  diffraction  pattern  with  a  large,  central  peak  known  as 
the  Airy  disk,  the  diffraction  radiation  pattern  produced 
by  electrons  is  characterized  by  a  null  at  0  =  0  and  is 
highly  dependent  upon  aperture  size.   In  addition,  a  strong 
peak  at  0  =  Y    is  present,  characteristic  of  transition 
radiation.   If  the  aperture  size  is  very  small,  the  structure 
in  the  diffraction  region  disappears  and  the  radiation 
pattern  approaches  that  of  pure  transition  radiation.   The 
only  radiation  present  is  the  strong  peak  at  0  =  y   •   I^~ 
creasing  the  aperture  size,  or  more  properly,  the  dimensionless 
parameter  k.n.    ,  decreases  the  peak  at  0  =  y    and  compresses 
the  structure  in  the  diffraction  region.   Eventually,  as  kn.^ 
approaches  infinity,  the  radiation  goes  to  zero.   The  diffrac- 
tion radiation  pattern  in  the  diffraction  region  exhibits 
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some  dependence  on  electron  energy  but  not  nearly  as  strong 
as  in  the  transition  region.   The  radiation  intensity  in 
both  the  transition  and  diffraction  regions  has  a  minor 
dependence  on  bunch  size. 

Ter-Mikaelian  describes  diffraction  radiation  as  having 
a  strong  peak  at  6  =  y    followed  by  a  broad  series  of  peaks 
and  nulls  dependent  on  J   (k  sin9l  ) .   The  diffraction  pat- 
terns developed  in  this  paper  are  primarily  dependent  on 
Jn  (k  sin0^  )  and  J„  (k  sin  dn.  ),   but  when  plotted  are  found 
to  be  in  complete  agreement  with  Ter-Mikaelian ' s  results. 
Both  results  are  in  agreement  on  several  key  aspects  of 
diffraction  radiation: 

1)  There  is  a  strong  dependence  of  diffraction  radiation 
on  aperture  size. 

2)  The  diffraction  radiation  energy  goes  to  zero  as 
aperture  size  goes  to  infinity. 

3)  The  diffraction  radiation  energy  is  dependent  on  the 
energy  of  the  beam,  described  by  y. 

4)  The  diffraction  radiation  exhibits  a  strong  peak  at 

6  =  y-l. 
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APPENDIX  A 

EVALUATION  OF  LINE  CHARGE  APPROXIMATION 
FOR  CYLINDRICAL  BUNCFI 

As  described  in  Chapter  II,  the  static  fields  for  a 
cylindrical  charge  distribution  were  approximated  by  a 
modified  line  charge.   The  linear  charge  distribution  was 
given  by 


X   =   PTT^Q  C/cm  , 


where 


a    =   cylindrical  bunch  radius  in  cm 
o       ^ 

p    =   charge  density'  of  cylindrical  charge 
bunch  in  C/cm-^  . 


When  b,  the  off-axis  distance,  is  less  than  a   the  linear 
'  o 

2 

charge   density   was    given   by    A    =    pTrb    . 


Y'     ^L 


pa.2i 


For  the  geometry  shown  above,  the  field  for  the  line 
charge  is  given  by 
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P"^o  1  ,     Z  Z-L 


4"^  b  >,^2  ^^2,1/2    ,^2  ^,2_^,2jl/2 


As  a  check  of  this  approximation,  the  fields  predicted 
by  the  linear  approximation  were  compared  to  the  fields 
resulting  from  the  actual  cylinder  of  charge,  which  were 
computed  by  numerically  integrating 


r     o            2tt           L                        X  -  r  cos  6  ^    ^a  ^r, 

j  (  (       5 5 ^-^^Ty  r  drde  dZ 


J  J  r  ^Y--r  r-ricfil     +^T-   c-infl^        +     (7.-7.' )^  \^/'^ 


47Te  ^Q         e=0         z=0  ((X-rcosG)   +(rsine)      +    (Z-Z')    ) 


Figure  A. 1.1  shows  the  comparison  for  three  "long" 
cylindeis  (L  >>  a  ) ,  where  the  observation  point  is  outside 
of  the  cylinder  (b  >  a  ) .   The  linear  approximation  is 
expected  to  work  best  in  this  situation. 

As  can  be  seen  from  the  figure,  the  approximation  is 
very  close  for  the  cylinders  of  lengths  1  cm  and  2  cm  and  is 
identical  for  the  longer  10  cm  bunch.   The  radial  fields 
become  very  nearly  equal  to  that  of  an  infinitely  long  line 
of  charge  very  close  to  the  ends  of  the  bunch. 

The  linear  approximation  is  expected  to  perform  worst 
for  a  relatively  short,  thick  cylinder  of  charge  when  the 
observation  point  is  on  the  interior  of  the  cylinder. 

Figure  Al . 2  shows  the  comparison  for  these  cylinders 
ranging  in  relative  length  from  L  =  a   to  L  >>  a  .   For  the 
shorter  length,  the  approximate  fields  are  found  to  rise 
slightly  more  sharply  and  to  higher  levels  than  the  actual 
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Figure    8:       E    (''"'^'"'/'^  )    vs    Z    for    cylinders    of    various    lengths 


a^  =    1    cm,    b=     .5    cm 
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fields.   Even  in  this  most  extreme  example,  however,  the 
approximation  should  be  adequate  for  our  purposes.   As  the 
bunch  length  increases,  the  actual  and  approximate  fields 
are  found  to  agree  very  closely,  with  a  slightly  steeper 
rise  in  the  linear  approximation. 

In  conclusion,  we  find  that  treating  the  cylindrical 
charge  distribution  as  a  modified  line  charge  offers  a  very 
satisfactory  approximation  to  the  actual  fields  present. 
Using  this  approximation  for  a  static  charge  bunch,  one  can 
compute  very  closely  the  Fourier  components  of  a  relativistic 
charge  bunch. 
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APPENDIX  B 
DIFFERENCES  WITH  THE  THEORY  OF  TER-MIKAELIAN 

In  Chapter  III,  it  was  brought  to  the  reader's  attention 
that  the  equations  derived  in  this  paper  were  slightly  at 
odds  with  the  results  presented  by  Ter-Mikaelian  in  Reference 
1.   In  his  paper,  Ter-Mikaelian  stated  his  results  for  an 
electron  passing  through  the  center  of  a  circular  aperture 
of  radius  a   for  ooa  /vy  <<  1  (in  Gaussian  units)  . 


t:;—      o     — T  J    i'^^      cos    \b 
2         Z    ,     Z      o     -    o 


2Tr    c    q     +a 


^ — -  J     (q/L    )  sin    ip 


\  ^2         2,20-^0 

-*  ZTT    c    q      +  a 


where : 


a      =      (jo/vy 

q      =      k  sine      =      projection   of   the  wave   vector   on 

the    z    =    0    plane 

\p      =   the  angle  between  q  and  the  X-axis 
These  equations  give  a  radial  field  equal  to 

E,   =  -^ 0    ^     n    J^(q-^^)   • 

-^      27T  c  q   +Y 

Since  this  equation  appears  to  be  inconsistent  with  the  | 

equations  derived  in  this  paper,  this  appendix  will  explain 
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at  length  the  method  used  to  obtain  the  Bessel  function  in 
the  equation. 

The  integral  over  d6  in  the  Fraunhofer  diffraction  problem 
was  written  as 


2tt 
f  ,       1  (kbX/R)cos  e  .„  ,, 

j         cos  9  e      •^  '      d0  db 


0 


The  cos  B  term  in  this  integral  is  due  to  the  radial  nature 

i  ft 
of  the  electric  field.   Using  the  expression  cos  0  =  (e    + 

—1  fl 
e   )/2,  the  integral  can  be  re-written  as 


2u   e^9  ^^    ^^   i(kbX/R)cos  9 ^^ 
0  ^ 


2tt   i(e+^^^os9)     i(-9+^^^os0) 
=   y   J    e  +  e  d( 

0 


.  -m    27T  .  ,   ,      . 
-,    ,    .         1     f         1  (mv+ucosv)  -, 
Since  J  (u)  =  -^ —     e  dv,  or 

m  ^'^    C] 


„     .m^     ,    ,  (■  i  (mv+ucosv)    dv    , 

2tti   J    (u)       =  e 

m  Q^ 


the  integral  over  d9  is  equivalent  to 


—1  k yh         kxb 

TTiJ^(kXb/R)  +  TTi  -"j  (kXb/R)   =   ui(j-^(^)  -  J_i(-p-) 
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Since  J'(x)  =  {j  _-,  (x)  -  J  _^-,(x)}/2,  with  n  =  0,  this 
expression  reduces  to 


-2.i  J'(S|b, 


which  equals 


2      , kXd , 
7T1  J^  (-R-) 


since  J'  (x)  =  -J-,  (x)  . 

As  shown  earlier,  the  sin  6  dependence  of  the  diffraction 
pattern  is  given  by 


2,zo-l    o2^o       l^o2    o 

q   +a 


This  appears  to  be  quite  different  than  that  found  by  Ter- 
Mikaelian . 

However,  if  one  plots  the  two  solutions,  he  will  find 
they  are  identical. 
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